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Thermodynamic properties of the quantum Heisenberg spin chains with S = 1/2, 1, and 3/2 are 
investigated using the transfer-matrix renormalization-group method. The temperature dependence 
of the magnetization, susceptibility, specific heat, spin-spin correlation length, and several other 
physical quantities in a zero or finite applied field are calculated and compared. Our data agree well 
with the Bethe ansatz, exact diagonalization, and quantum Monte Carlo results and provide further 
insight into the quantum effects in the antiferromagnetic Heisenberg spin chains. 
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I. INTRODUCTION 

There are a lot of quasi-one dimensional compounds 
whose behaviors can be adequately described within the 
framework of interacting spin chains governed by the 
Heisenberg model. Extensive studies on this model have 
shed light on the quantum nature of spin dynamics. 
In 1983, Haldane predicted that the one-dimensional 
Heisenberg antiferromagnetic model with integer spin has 
an excitation gap and a finite correlation length [jlj . Since 
then a great amount of experimental and theoretical ef- 
fort has been made towards understanding the difference 
between half- integer and integer spin chains. 

In this paper we report results of a transfer-matrix 
renormalization-group (TMRG) study for the ther- 
modynamics of quantum Heisenberg spin chains. Wc 
have calculated the magnetic susceptibility, specific heat, 
spin-spin correlation length, and other experimentally 
relevant quantities as functions of temperature and ap- 
plied magnetic field for the S=l/2, 1, and 3/2 spin chains. 

The Heisenberg model is defined by the Hamiltonian 



H 



N 
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where J is the spin exchange constant and H is an ap- 
plied magnetic field In this paper, we consider an- 
tiferromagnetic spin chains only. We use units in which 
J = 1. 

The spin-1/2 Heisenberg model is integrable by Bethe 
ansatz. Many of its thermodynamic quantities, for exam- 
ple, the specific heat and the spin susceptibility, can be 
calculated by solving the Bethe ansatz equations. The 
conformal field theory is also very useful in analysing 
the low temperature low field thermodynamic properties, 
since the S=l/2 Heisenberg model is equivalent to the 
k=l Wess-Zumino-Witten nonlinear a model. 

Higher spin Heisenberg models are not at present 
amenable to rigorous solution. The finite temperature 
properties of the model were studied mainly through 
transfer matrix ||, quantum Monte Carlo and 
other approximate method jlO 11|. 



At zero field, the ground state of H is a spin sin- 
glet with zero spin magnetization. At finite field, the 
magnetization of the ground state becomes finite and 
increases with increasing H. There is a critical field 
H C 2 beyond which all spins are fully polarized at zero 
temperature. If we denote E(N, S) as the lowest eigen- 
value of a N-site Heisenberg chain with total spin S at 
H = 0, then it is straightforward to show that H C 2 = 
E(N, S max ) - E(N, S max - 1) = AS, independent on N. 
Below H C 2, a canted Neel order, namely a state which 
has both ferromagnetic order along the z-axis and anti- 
ferromagnetic order in the xy plane, exists at sufficiently 
low temperature, and the pitch vector (i.e. the value of 
the momenta at which the static structure factor shows a 
peak) decreases continuously from ir to with increasing 
H. 

Integer spin chains can be described by the quantum 
nonlinear a model. It is from the study of this model 
that Haldane made the famous 'Haldane conjecture'. The 
ground state of the 0(3) a model has a finite excitation 
gap and consequently a finite correlation length. The ap- 
plication of a magnetic field causes a Zeeman splitting of 
the triplet with one member crossing the ground state at 
a critical field H c \ whose value is equal to the excitation 
gap A. When H < H cl , the Haldane gap persists and the 
ground state is still a non-degenerate spin singlet state. 
When H c i < H < H C 2, the ground state has a nonzero 
magnetization with gapless excitations. Thus across H c i, 
an integer spin system undergoes a commensurate to in- 
commensurate transition. This is an interesting feature 
which is absent in a half odd integer spin system. Evi- 
dence for such transitions has been used to identify the 
Haldane gap in NENP |l2j and other quasi- Id spin com- 
pounds 1 1 3 1 - Just above H&, the ground state can be 
regarded as a Bose condensate of the low energy boson. 
Varying the magnetic field is equivalent to varying the 
chemical potential for this boson and the (uniform) mag- 
netization corresponds to the boson number JlJ] . 
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II. TMRG 

In this section, we discuss briefly the TMRG method. 
A more detailed introduction to the method can be found 
from Refs. |,|. 

The TMRG method starts by mapping a Id quantum 
system onto a 2d classical one with the Trotter-Suzuki 
decomposition and represents the partition function as a 
trace of a power function of virtual transfer matrix Tm'- 
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where M is the Trotter number. Tm is defined by an 
inner product of 2M local transfer matrices 
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and t — (3/M. \af) is an eigenstate of 5f and erf is 
the corresponding eigenvalue: $f|of ) = of The su- 
perscripts and subscripts in Tm and £ represent the spin 
positions in the Trotter and real space, respectively. 

Tm conserves the total spin in the Trotter space, i.e. 
Efc °f • Thus Tm is block diagonal according to the value 



of Ea 



For the S=l/2 Heisenberg model, it was 



shown rigorously that the maximum eigenstate of Tm is 
non-degenerate and in the J2k a i = subspace, irrespec- 
tive of the sign of J and the value of H [[15| . When S > 
1/2, we found numerically that the maximum eigenvec- 
tors of Tm are also in the J2k a i = subblock. 

In the thermodynamic limit, the free energy per spin, 
is given by 
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lim 
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where X ma x is the maximum eigenvalue of Tm- From the 
derivatives of F one can in principle calculate all ther- 
modynamic quantities. The internal energy U and the 
spin magnetization M z could, for example, be calculated 
from the first derivative of F with respect to H and T, 
respectively. However, numerically it is better to calcu- 
late U and M z directly from the eigenvectors of Tm [@ • 
The spin susceptibility \ — dM z /dH and the specific 
heat C = OU/dT can then be calculated by numerical 
derivatives. The specific heat such determined is gen- 
erally found to be less accurate than, for example, the 
susceptibility data at low T. The reason for this is that 
U changes very slowly with T (or equivalently C is very 
small) at low T, and a small error in U would lead to a 
relative large error in C. 

The correlation length of the spin-spin correlation 
functions, defined by = — limL_ >00 ln(5i. ct S , i+L,a), can 



also be calculated from this method. The longitudinal 
and transverse correlation lengthes are determined by 
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where A2 is the second largest eigenvalue of Tm in the 
subspace J2k a l = an< ^ ^1 is the largest eigenvalue of 



(2) Tm in the subspace J2k a l = ^1- 



Figure [l] shows the configuration of superblock used 
in our calculation. This configuration of superblock is 
different from those used in Refs. I^Q. The advantage for 
forming the superblock in such a way is that the transfer 
matrix Tm in this case can always be factorized as a 
product of two sparse matrices (which are block diagonal 
with respect to n s and as <8 n e , respectively). To treat 
these two sparse matrices instead of Tm itself allows us 
to save both computer memory space and CPU time. 
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n P 



FIG. 1. A pictorial representation of a superblock, which 
consists of system and environment blocks connected by three 
added spins to form a periodic chain. n s and n e are indices for 
the system and environment blocks, respectively. Initially the 
system block contains two spins and the environment block 
contains only one spin. At each iteration the system and 
environment blocks are augmented by adding a 1 to each of 
the blocks. The augmented system h s is formed by n e a 1 , 
and the augmented environment n e is formed by a 1 ® n s . 

We compute the maximum eigenvalue, \ m ax, and the 
corresponding right and left eigenvectors, \ip R ) and (ip \, 
of Tm using an implicitly restarted Arnoldi method [ 16) . 
This method is more efficient than the power method 
which we used before HQ- Ai,2 can also be calculated 
from this method, but their truncation errors are gen- 
erally larger than those of X ma x- Thus the correlation 
lengthes determined from Eqs. (|^, ||) are generally ex- 
pected to be less accurate than the free energy or other 
thermodynamic quantities. 

In the TMRG method, the density matrix for the aug- 
mented system or environment block is non-symmetric, 
which is different than in the zero-temperature DMRG 
method. Numerically it is much more difficult to treat 
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accurately a nonsymmetric matrix than a symmetric one 
because the errors in \4j r ) and (fp L \ may affect the the 
(semi-) positiveness of the density matrix and increase 
the truncation error of the TMRG. 

The TMRG treats directly an infinity spin chain. 
There is therefore no finite size effect. The error caused 
by the finiteness of the Trotter number (or r) is of order 
t 2 , which is generally very small. The error resulted from 
the truncation of basis states is difficult to estimate. A 
rough estimation for this type of error can be obtained 
from the value of truncation error, which is smaller than 
10~ 3 in all our calculations. More accurate results can be 
obtained simply by extrapolating the results with respect 
to both t and the number of states retained m. 



III. RESULTS 
A. S=l/2 

The S=l/2 Heisenberg model is by far the best under- 
stood spin system. In the absence of field, the ground 
state is massless and the Bethe ansatz result for the 
ground state energy is (1/4 — In 2). The lowest exci- 
tatins states are spin triplets with a spin wave spectrum 
jlTj e(k) — (7r/2) I sinfc|. Above this lower boundary of 
excitations, there is a two-parameter continuum of spin 
wave excitations with an upper boundary given by |l§| ] 
e(k) = 7r| sin(/c/2)|. There are other excitations above 
this upper boundary. 

The specific heat of the model was first calculated nu- 
merically by Bonner and Fisher |l9|. They found that 
at low temperature C ~ 0.7T when H = 0. Later, Af- 
fleck ^0|, using the conformal field theory, found that 
C = (2/3)T. The zero field zero temperature magnetic 
susceptibility is x = 1/tt 2 - This result was first obtained 
by Griffiths with the Bethe ansatz and numerics | pl| . 
Recently, Eggert et al [^2| found that the zero field sus- 
ceptibility approaches to its zero temperature value log- 
arithmically when T < 0.01, 
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where v = (tt/2) is the spin wave velocity and Tq « 7.7. 

Figure |] shows the TMRG results for a number of ther- 
modynamic quantities of the S=l/2 antiferromagnetic 
Heisenberg model in various magnetic fields with to = 81 
and r = 0.1. At zero field, the TMRG reproduces accu- 
rately the results which were previously obtained by the 
Bethe ansatz or conformal field theory. The extrapolated 
zero field zero temperature values of the internal energy 
U(T) (i.e. ground state energy), the spin susceptibility 
x(T) and the linear coefficient of the specific heat C(T), 
are -0.443, 0.109 and 0.66, respectively. In all the fields 
which we studied, the peak position of x{T) is located at 



a temperature which is about twice of the peak temper- 
ature of C{T)/T. This is due to the fact that x(T) is a 
measure of two-particle excitations and C(T)/T is only 
a measure of one-particle density of states. The maxi- 
mum of x when H = is approximately equal to 0.147 
at T = 0.64, consistent with the Bethe ansatz calculation 

©■ 

There is a significant change in the temperature depen- 
dences of x an d C when H is below and above a critical 
field H c2 = 2. Below H c2 , both C/T and x are finite at 
zero temperature, which shows that gapless spin excita- 
tions with a finite low energy density of states exist in this 
regime. When H = H C 2, both x and C/T vary as 1/vT 
at low temperature; the extrapolated value of VT~X an d 
C/VT at zero temperature are 0.152 and 0.22, respec- 
tively. The divergency of x an d C/T at T = implies 
that the density of states of spin excitations is divergent 
at zero energy when H = H C 2- Above H C 2, there is a 
gap in the excitation spectrum as both x & n d C/T drop 
to zero exponentially at low temperatures. The value of 
the gap estimated from the low temperature behavior of 
X and C grows linearly with H — H c2 . 
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FIG. 2. (a) Free energy F(T), (b) internal energy U(T), 

(c) ratio between the magnetization M Z (T) and the field H, 

(d) spin susceptibility x(T), (e) specific heat C(T), and (f) 
C(T)/T, per lattice site, for the spin- 1/2 Heisenberg model 
in four applied fields, H = (or 0.05), 1, 2, and 2.5. When 
H = 2, T 1/2 x(T) and 2C(T)/T 1/2 , instead of X (T) and 
C(T)/T, are shown in (d) and (f), respectively, r = 0.1 and 
m — 81 are used in the TMRG calculations. 



Figure (||) shows the longitudinal and transverse cor- 
relation lengthes of the S=l/2 model in different applied 
fields. When H = 0, £ x — £ z diverges at zero temperature 
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and the slope of ^J 1 at low temperature is approximately 
equal to 2, in agreement with the thermal Bethe ansatz 
as well as the k=l WZW a model result [B3I 
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In the presence of magnetic field, £ z is substantially 
suppressed and becomes finite at zero temperature. As 
the correlation length is inversely proportional to the en- 
ergy gap of excitations, the finiteness of £ z at T = 
means that the longitudinal spin excitation modes are 
massive in a field. There is a small dip in the curve of 
£ z at T ~ 0.4 when H = 0.2 or at T ~ 0.6 when H = 1. 
This dip feature of £ z , as will be shown later, appears 
also in large S systems. 

The effect of the applied field on the transverse spin 
excitation modes is weaker than the longitudinal modes. 
The transverse spin excitations become massive only 
when H > H C 2- Below H C 2, £ x drops to zero linearly 
with T, as for the case H — 0. When H = H C 2, £ x 
drops to zero as VT. Clearly the thermodynamics of the 
Heisenberg model in a field is mainly determined by the 
transverse excitation modes at low temperatures. 




FIG. 3. (a) and (b) £ z vs T for the spin- 1/2 Heisen- 

berg models, r = 0.1 and m — 81 are used in the TMRG cal- 
culations. The circles and squares are thermal Bethe ansatz 
results. 



A simple understanding of the above results can be 
obtained from an equivalent spinless fermion model of 
the S=l/2 Heisenberg model: 

h = -k c Ui* + 4^+0 + (h- i)4* 



+(| - H) + c\ci& i+1 c i+1 
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where Cj a spinless fermion operator which is linked to 
the S = 1/2 spin operator by the Jordan- Wigner trans- 
formation Cj = Sf exp(«7r X)z<i c i The magnetic field 
is equivalent to a chemical potential for the fermions. 
When H=0, the ground state has zero uniform mag- 
netization, corresponding to a half filled fermion band. 
As H increases, the ground state is ferromagnetically 
polarized and the Fermi energy shifts down. When H 
is smaller than the critical field H c2 , the ground state 
of this fermion model has no gap but the spin ori- 
entation is canted, i.e. in an incommensurate state. 
The pitch angle of this incommensurate ground state, 
namely the wave vector at which the maximum of the 
spin-spin correlation function appears, can be estimated 
from the Fermi momentum of non-interacting fermions as 
2kF = [1 — 2M Z (H)]%. This value of &p is not normalized 
by interactions according to the Luttinger theorem. 

Above H C 2, there are no fermions at the ground state. 
At low temperature, the number of fermions excited from 
the ground state are rare. Thus, as a good approxima- 
tion, the interaction term in (||) can be ignored at low 
temperatures. For the non-interacting system, the en- 
ergy dispersion of fermion excitations from the ground 
state is given by ek — H — (1 + cos k), which has a gap of 
H — H c2 . When H = H c2 , £fc = (1 — cos k) , the density of 
states of excited fermions varies as l/y 7 ^ at low energy. 
From the standard theory of noninteracting fermions, it 
is straightforward to show that this singular density of 
states will cause both \ an d C/T to diverge as l/VT 
at low temperature. When H > H C 2, there is a gap in 
the fermion excitations, both \ an d C should decrease 
exponentially at low temperature. These results are just 
what we found in Figure 



B. S=l 

As mentioned previously, there are two critical fields in 
the spin-1 Heisenberg model, H c i and H C 2\ below H c \, 
the ground state is a massive spin singlet; above H c2 , 
the ground state becomes a fully polarized ferromagnetic 
state; between H c \ and H C 2, the ground state is massless 
and has a finite magnetization. The temperature depen- 
dence of thermodynamic quantities of the S=l model be- 
low, above and (approximately) at these critical fields is 
shown in Fig. |. 

At zero field both x(T) and C(T) drop exponentially 
with decreasing T at low temperatures duo to the open- 
ing of the Haldane gap. In this case the ground state 
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energy extrapolated from the internal energy U(T) is - 
1.4015, in agreement with the zero-temperature DMRG 
result . The ground state excitation gap A can be de- 
termined from the temperature dependence of \ and C at 
low temperatures. If we adapt the ansatz that the low- 
lying excitation spectrum has approximately the form 
M 



Fig. | shows y/2TC(T)/3 X (T) as a function of T for 
the S=l Heisenberg model at zero field. By extrapola- 
tion, we find that A ~ 0.41 in agreement with the zero 
temperature DMRG and exact diagonalization |2S|] 
results. Given A, the value of v can be found from Eq. 
( |Io| ) in the limit T — ► 0. The value of v we found is 
~ 2.45, consistent with other numerical calculations [E6| . 
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where v is the spin wave velocity, it is then straightfor- 
ward to show that, when T <gC A, the spin susceptibility 
and the specific heat are 
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irrespective of the statistics of the excitations. Taking 
the ratio between x(T) and C(T) gives 



A=limJ^2. 



(12) 



This is a very useful equation for determining A, espe- 
cially from the point of view of experiments, since both 
x(T) and C(T) are experimentally measurable quanti- 
ties. 
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FIG. 4. Thermodynamic quantities for the spin-1 Heisen- 
berg model in five applied fields, H = (or 0.05), 0.4105, 2, 
4, and 4.5. When H = and 0.4105, F(T) and £/(T) are 
almost indistinguishable in the figure. When H — 0.4105 or 
4, T 1/2 x(T) and 2G(T)/T 1/2 (dotted lines) are shown in (d) 
and (f), respectively, r = 0.1 and m = 100 are used in the 
TMRG calculations. 
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FIG. 5. (2TC(r)/3x(T)) 1/2 vs T for the S=l Heisenberg 
model, (m = 100 and r = 0.1.) 



At the two critical fields, H c i ~ 0.4105 and H C 2 = 4, 
both x(T) and C(T)/T diverge as T" 1 / 2 at low temper- 
atures. This divergence, as for the S=l/2 case, is due to 
the square-root divergence of the density of states of the 
low- lying excitations. At H c \, one branch of the S=l ex- 
citations becomes massless and the low-energy excitation 
spectrum is approximately given by 



e(k)~— (fc-Tr) 2 . 



(13) 



If we assume that these excitations are fermion-like, i.e. 
satisfy the Fermi statistics (a short-range interacting 
Bose system is equivalent to a system of free fermions), 
then it is simply to show that when T -C A 
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By extrapolation, the TMRG result gives 
0.45, which is close to that given by Eq. @. The value 
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of T~ 1 / 2 C\t^o is difficult to determine accurately from 
the TMRG result since the error of C(T) at low temper- 
ature is larger than C(T) itself. 
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FIG. 6. (a) and (b) CJ 1 vs T for the spin-1 Heisenberg 
models. ^(H = 0) and ^(H = 0.4105) are almost indis- 
tinguishable when T > 0.6. r = 0.1 and m = 100 are used in 
the TMRG calculations. 

Fig. H shows and ^J 1 of the S=l model in different 
applied fields. Because of the Haldane gap, the longitu- 
dinal correlation length is finite even at zero field. Thus 
the longitudinal excitation mode is always gapful. Above 
H c i but below H C 2, the transverse correlation length di- 
verges at zero temperature. Thus the transverse mode is 
massless in this field regime. We note, however, that 
at H — (similarly £~ 1 at H — 0.4105) does not decrease 
monotonically at low temperatures, actually it starts to 
raise below a temperature T* . This non-monotonic be- 
havior of at H = seems only due to the truncation 
error. Our preliminary calculation shows that T* de- 
creases with increasing m. But by far we still do not 
know if T* is zero in the limit m — > oo. Further inves- 
tigation to the low temperature behavior of £ is needed. 
If we do an extrapolation using the TMRG data of £~ 
above T*, we find that £, X (T = 0) ~ 6.0, which is con- 
sistent with the zero temperature DMRG result pjj. If, 
however, the TMRG data below T* are also included in 
the extrapolation, we find that f x (T = 0) is only ~ 4.4. 



When H > H c i, drops rather sharply at some 
temperatures. These sharp drops of £~ happen when the 
second and third eigenvalues of the transfer matrix with 
J2k <J k = cross each other. The physical consequence 
of these sudden changes in the longitudinal correlation 
length is still unknown. 

C. S=3/2 

The thermodynamic behaviors of the S=3/2 Heisen- 
berg model, as shown in Fig. [?], are similar to those of 
the S=l/2 model. When H < H c2 = 6, \ is always finite 
at zero temperature, indicating that the ground state is 
massless; above H C 2, the ground state is fully ferromag- 
netic polarized and a gap is open in the excitation spec- 
trum; at H c2 , both x(T) and C/T diverge as T^ 1 / 2 . At 
zero field, the susceptibility data, extrapolated to zero 
temperature, gives x(0) ~ 0.67, consistent with recent 
numerical calculations |pL|ll). 
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FIG. 7. Thermodynamic quantities for the spin-3/2 
Heisenberg model in four applied fields, H = (or 0.05), 3, 
6, and 6.5. When H = 6, T 1/2 X (T) and 2C(T)/T 1/2 (dotted 
lines) are shown in (d) and (f), respectively, r = 0.075 and 
m = 81 are used in the TMRG calculations. 

The crossover from quantum to classical behavior can 
be clearly seen (Fig. |^) by comparing the zero field sus- 
ceptibility and specific heat of the S=l/2, 1 and 3/2 
spin chains with the corresponding results of the clas- 
sical Heisenberg spin chain: J32| 

,_. 1 1 — u(T) ,. 1 . . 

X(T) = 7=7. u(T)=coth T, (18) 

AV ; 3Tl + u(T)' v ; T y ' 
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C(T) = 1 - 
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T 2 smh z (l/T) 



(19) 



At high temperatures, T/S(S + 1) > 1, the quantum re- 
sults approach asymptotically to the classical ones. The 
agreement between the quantum and classical results per- 
sist down to progressively lower temperatures as S in- 
creases. At low temperatures, however, the difference 
between the results of the S=3/2 system and those of 
the classical model is still very large, indicating the im- 
portance of the quantum effects in the study of quantum 
spin chains. (Note for the classical Heisenberg model, 
C(T) docs not vanish at zero temperature. This is a 
unrealistic feature of this model.) 
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8. The zero field specific heat and susceptibil- 
ity vs T/|S| 2 for the S=l/2, 1, and 3/2 spin chains 
(|S| 2 = S(S + 1)). The corresponding results for the classical 
Heisenberg model (|S| 2 = 1) are also shown for comparison. 



D. Zero temperature magnetization 

Figure || shows M Z (T — 0), extrapolated from the fi- 
nite temperature TMRG data of M Z (T), for the S=l/2, 
1, and 3/2 systems. For comparison the Bethe ansatz 



result |2lJ] for the S=l/2 Heisenberg model is also shown 
in the figure. With increasing S, we found that M z tends 
to approach to its classical limit (S — > oo) where M z 
increases linearly with H, i.e. M z /S = H/H C 2. 
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FIG. 9. Normalized zero temperature magnetization M z 
as a function of H/H c2 . The TMRG results for S=l/2 
(circles), S=l (solid line), and S=3/2 (dashed line) are ob- 
tained by extrapolation from the low T values of M z with 
(m,r) = (81,0.1), (81, 0.1), and (60, 0.75), respectively. Dot- 
ted line is the Bethe ansatz result for the S=l/2 system. 



For the S=l/2 system, the TMRG result agrees well 
with the Bethe ansatz one. A least square fit to the 
curve of M z up to the fourth order term of \/H C 2 — H 
gives M = M c — a lv ?H c2 -H + a 2 (H c2 -H)- a 3 {H c2 - 
Hf' 2 + ai(H c2 - H) 2 , with ax w 0.448, a 2 w 0.123, 
as = 0.05 and — 0.00744. a\ agrees very accurately 



with the exact value V2/tt. In the weak field limit, 
asymptotic behavior of M z is 



the 



M z « Z ( I 



2ln(H/H c2 ) 



(20) 



as predicted by the Bethe Ansatz theory |29|j30|] . 

For the S=l model, M Z {T — 0) becomes finite when 
H > H c \. In a very narrow regime of field near H c \, 
M z {0) varies approximately as y/H — H c \ , in agreement 
with the prediction of the Bose condensation theory jidj . 
But the difference between the result of the Bose conden- 
sation theory 



M Z {T = 0) 



sj2{H - H cl )A 



TV I.' 



(21) 



and that of the TMRG becomes already significant when 
H - H cl = 0.04. 

Near H c2 , M z approaches to its saturation value M c = 
S as a function of y/H c2 — H for the three spin systems 
we study. For the S=l/2 and 1 systems, the TMRG 
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results agree very accurately with the Bethe ansatz result 

MM 



S 



2 



H 

H C 2 



(22) 



For the S=3/2 system, the asymptotic regime of H is 
very narrow, we cannot do a detailed comparison between 
the TMRG result and Eq. (p2|). The magnetization curve 
does not show a plateaus at M z = 1/2, in agreement with 
other studies fl33l]. 



according to a scaling argument [B4| . Our TMRG result 
agrees well with this equation at low temperatures. By 
plotting (Txs) 2 against InT, we find that D x w 0.30 and 
Ty w 10.5, in agreement with a Monte Carlo calculation 



For the S=l model, Xs is finite at zero temperature. 
The extrapolated zero temperature value of Xs for the 
S=l model is 18.55, in agreement with other numerical 
calculations fi^Hl. 



IV. CONCLUSION 



E. Staggered susceptibility 

To calculate the staggered susceptibility, we add a stag- 
gered magnetic field H s to the Hamiltonian H . The stag- 
gered magnetization is then calculated in a way similarly 
to the calculation of the uniform magnetization. The 
staggered susceptibility %s is obtained by differentiating 
the staggered magnetization with respect to H s . For half- 
odd-integer spin chains Xs diverges as T~ l at low temper- 
atures. Thus the staggered magnetization becomes satu- 
rated at low temperatures when Xs(T)H s > S, since the 
maximum value of the staggered magnetization per spin 
is S. This means that to estimate accurately the zero field 
staggered susceptibility the staggered field used should 



satisfy the condition H s <C S/xs{T„ 
the lowest temperature to study. 



,), where T„ 



S=3/2 




0.1 h , L 

0.01 0.1 1 10 

T/S(S+1) 

FIG. 10. the staggered susceptibility Xs(T) at zero field. 
The parameters used for the S=l/2, 1, and 3/2 chains are (H s , 
t, m) = (0.0001, 0.1, 140), (0.0001, 0.05, 100), and (0.0001, 
0.025, 81), respectively. 



In conclusion, the temperature dependence of the sus- 
ceptibility, specific heat, and several other quantities of 
the quantum Heisenberg spin chains with spin ranging 
from 1/2 to 3/2 in a finite or zero applied magnetic field 
are studied using the TMRG method. At high temper- 
atures, the quantum results for the specific heat as well 
as other thermodynamic quantities approach asymptoti- 
cally to the classical ones for both the integer and half- 
integer spin systems. At low temperatures, however, the 
quantum effect is strong and the integer spin chains be- 
have very differently than the half- integer spin chains. 
For the S=l model, both x an d C decay exponentially 
at low temperatures due to the opening the Haldane gap. 
For the S=l/2 and 3/2 spin chains, there is no gap in the 
excitation spectrum and both x an( i C/T are finite at 
zero temperature. The thermodynamics of the Heisen- 
berg spin chains in an applied field is mainly determined 
by the transverse excitation modes. At low temperatures, 
X, C/T, and the transverse correlation length £ x diverge 
as T- 1 / 2 at H c2 for the S=l/2 and 3/2 models and at 
both H c i and H C 2 for the S=l model. This square-root 
divergence indicates that the low energy spin excitations 
have a square-root divergent density of states at these 
critical fields. Our data agree well with the Bethe ansatz, 
quantum Monte Carlo, and other analytic or numerical 
results. 
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Fig. [To] shows the temperature dependence of Xs for 
the S=l/2, 1, and 3/2 spin chains. At low temperatures, 
Xs for the S=l/2 model is expected to have the following 
asymptotic form 

Xs = ^^fJf (23) 
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